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FOREWORD 


This report is one of several to be published from research 
conducted under, or supported in part by, NASA Contract’ No. NAS8- 
26751, entitled "Atmospheric Energetics as Related to Cyclogenesis 
Over the Eastern United States." A number of approaches have been 
and continue to be followed in the conduct of the research. The 
results presented in this report are those from one approach, and 
represent only a portion of the total research effort. Other 
reports will be published as the research progresses. 
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INTRODUCTION 


a. Statement of the problem 

While meteorologists have long realized that the migratory 
extratropical disturbances are areas where energy conversions are 
maximized, the exact mechanisms that effect the conversions between 
kinetic, internal, and gravitational potential energy have remained 
somewhat obscure. In particular, the determination of how and 
where these conversions take place in a developing cyclone is still 
a fundamental problem in meteorology. 

Comprehensive studies of the energy processes in individual 
synoptic systems have been limited primarily because of the dif- 
ficulties in evaluating the boundary integrals that can represent 
a significant portion of the energy budget. The existence of non- 
zero boundary integrals also complicates the interpretation of the 
energy conversion processes. These factors have resulted in a 
majority of the effort being directed toward specifying the global 
energy budget, a simpler problem because the boundary integrals 
become zero when the entire atmosphere is under consideration. 

Quantitative studies of the energy transformations in indi- 
vidual cyclones are necessary in order to better understand their 
evolution. For example, it is desirable to learn more about the 
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energy transformations necessary to maintain the kinetic energy 
of the cyclone against the ever present frictional dissipation. 
Admittedly, energy considerations represent only one aspect of the 
problem, and a complete understanding of the initiation, growth, 
and subsequent decay of extratropical cyclones would involve the 
simultaneous solution of the mass, momentum, and energy relation- 
ships. However, this is an unsolvable problem with the existing 
knowledge and observational network. 

The hypothesis underlying the research presented in this 
paper is that energy processes can be related to cyclogenesis 
through a detailed study of the energetics of a limited atmospher- 
ic volume that encompasses the cyclone. The validity of this 
hypothesis depends upon the extent to which: 

1) Energy processes can be defined in a limited region of the 
atmpsphere, 

2) The variability of the energy processes can be specified 
in both time and space. 

3) Boundary fluxes of energy can be determined. 

4) The interaction of the cyclone with its environment can 
be determined. 

The method that was used to study the energy budget of extra- 
tropical disturbances in the present research represents an 
extension of previous investigations in several respects. First, 
all the terms of the energy budget equations were considered, 
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including boundary fluxes of energy that were ignored in most 
previous studies. Second, objectively analyzed data were used in 
the evaluation of the various terms of the energy budget equations. 
Third, the energy budget equations were evaluated throughout the 
spatial and temporal extent of the cyclone. Results of previous 
studies were presented in terms of averages, thus providing no 
information on the variation of the energy processes in time and 
space. Finally, the interaction of the cyclone with the larger 
scales of motion is evaluated in a unique manner. While approaches 
other than the one presented in this paper certainly are feasible, 
this approach has the advantage of describing the energy budget in 
terms of the basic physical processes responsible for the energy 
transformations . 

b. Related studies 

The first description of a possible means of generating atmo- 
spheric motions is credited to Margules (Hess, 1959). He considered 
two adjacent air masses of different density that were allowed to 
adjust to a statically stable density stratification which was 
invariant in the horizontal direction - often referred to as the 
reference state. During this readjustment, potential energy is 
converted into kinetic energy. Margules' computations showed that 
the wind speed generated by this process is comparable to that 
observed in the atmosphere. 
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Margules' results, however, apply only to closed atmospheric 
systems, i.e., systems that include the entire atmosphere. Spar 
(1950), in a study of energy conversions in a developing cyclone, 
demonstrated that increases in the kinetic energy of open systems 
are not necessarily accompanied by decreases in potential energy, 
as Margules' theory predicts. This is because the boundary fluxes 
of energy are important in open systems and must be included in the 
energy budget. An extensive treatment of the role of boundary 
conditions in open systems has been given by Miller (1950; 1951) 
and Van Mieghem (1951) . It has become clear, however, that the 
rising of warm air accompanied by the sinking of cold air, as 
described by Margules, is one of the most important mechanisms 
responsible for the generation of kinetic energy in a developing 
cyclone. Margules also initiated the concept of available poten- 
tial energy, which is defined to be the difference between existing 
total potential energy (internal + gravitational potential) and 
the total potential energy that would exist after the entire atmo- 
sphere undergoes an isentropic redistribution of mass to the 
reference state. 

With the definition supplied by Margules, Lorentz (1955) 
developed an approximate expression for the generation of available 
potential energy. The definition of available potential energy, 
however, is based on the redistribution of mass of the entire atmo- 
sphere toward the reference state, and, until recently, 
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investigators have found it impossible to apply rigorously the 
available potential energy concept to a limited region of the 
atmosphere. This task has been accomplished by Johnson (1970) 
who has succeeded in prescribing an exact mathematical theory for 
available potential energy of open systems. Diagnostic studies 
of energy conversions in developing cyclones based on the theory 
provided by Johnson should prove helpful in understanding the 
relationships between energy conversions and cyclogenesis. The 
concept of available potential energy is useful because it isolates 
that portion of the total potential energy which is available for 
conversion to kinetic energy. On the other hand, there is no 
assurance that all of this energy will in fact be converted, and 
some controversy exists over what constitutes a realistic state. 

Two additional sources for the supply of kinetic energy 
found in extratropical disturbances' were presented in classic 
papers by Charney (1947) and Kuo (1949) . Charney investigated a 
baroclinic atmosphere and found that the unstable growth of the 
perturbation in a baroclinic atmosphere results from a conversion 
of potential energy into kinetic energy. The instability is 
related to vertical wind shear of the basic current, temperature 
lapse rate, latitude, and wave length. In an extension of this 
work, Thompson (1961) has shown that the increase in the kinetic 
energy of the perturbation is maximized when the thermal wave 
lags the pressure wave by 90 deg. 
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Kuo investigated a barotropic atmosphere and found that 
perturbations in the basic flow will grow when a critical horizon- 
tal shear exists. The kinetic energy of the perturbation is 
increased at the expense of kinetic energy of the basic current. 
This process represents a transfer of existing kinetic energy from 
one state to another rather than a conversion of energy from one 
form to another. 

Although the theories of Kuo and Charney appear conflicting, 
in reality both mechanisms are important and probably occur simul- 
taneously during cyclone development. The relative importance of 
the two mechanisms is uncertain; although, in a study of one 
particular cyclone (Johnson, 1970), it was found that two- thirds 
of the kinetic energy was supplied by the conversions from poten- 
tial energy and one-third by the kinetic energy of the basic flow. 

Since the appearance of these classic papers, innumerable 
articles treating various aspects of energy transformations have 
appeared in the literature. These studies can be divided conve- 
niently into two categories. The first category is concerned with 
average conditions over a long period of time and over a large 
portion of Earth's surface. This type of study is useful for a 
better understanding of the general circulation. A summary of a 
significant number of these articles and the results has been 
presented by Oort (1964) . 
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The second class of studies deals with kinetic energy gener- 
ation within individual synoptic systems. The work by Spar (1950) 
that has already been mentioned is one example. Other notable 
examples are Palmen (1958), Palmen and Holopainen (1962), Danard 
(1964), and Eddy (1965). In all of these studies, except for the 
one by Spar, the term has been used to represent the con- 

version rate between available potential energy and kinetic energy. 
In the above integral, go is the vertical velocity in pressure 
coordinates, o l is specific volume, V is the volume over which the 
integration is performed, and v" is a dummy variable of integration 
corresponding to V. When warm air is rising and cold air is sinking, 
the term will be positive, indicating an increase in kinetic 

energy (White and Saltzman, 1956). However, because of the pres- 
ence of boundary fluxes of energy, there is no assurance that this 
term represents the correct value for local conversion rates in 
open systems (Sechrist and Dutton, 1970; Smith and Horn, 1969). 

This point will be clarified in a later discussion. 

Through the use of the cOit-technique, Palmen (1958) investi- 
gated the development of hurricane Hazel into an extratropical 
cyclone, and found that the conversion rate amounted to 51 watts 
m during the period of strongest development. In similar studies, 

but for separate synoptic cases, Palmen and Holopainen (1962) 

_2 

computed the conversion rate to be 21 watts m ; Danard (1964) 

_ o 

found a value of 16 watts m ; and Eddy (1965) found that the most 
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_ o 

intense conversion rate, 10 watts m , occurred in the middle 
troposphere ahead of the storm and during the initial development. 

If the case involving hurricane Hazel is disregarded, the 
conversion rates described above are in general agreement. The 
development of hurricane Hazel into an extratropical cyclone was 
an unusual event, and cannot be compared with typical extratropical 
cyclone development. 

A unique approach to the study of energy conversions in 
individual cyclones was developed by Sechrist and Dutton (1970) . 
They considered volumes that moved isentropically with the wind, 
and thus were able to eliminate boundary conditions and make more 
precise measurements of the rates of kinetic energy conversion 
within the volume. In addition, they were able to deduce the 
geographical sources of the kinetic energy, something not previous- 
ly possible. Since their results were based on the study of one 
cyclone, no general conclusions can be made. This approach is 
promising, but is applicable only when diabatic effects are at a 
minimum. 

After the development of the baroclinic instability theory by 
Charney, the popular opinion evolved that the extratropical cyclone 
was maintained against frictional dissipation by the release of 
this baroclinic instability, and diabatic effects were considered 
to be minor. Petterssen (1956), however, stressed the physical 
significance of geographical heat sources and sinks in intensifying 
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cyclonic or anticyclonic circulations. Further work by Palmen 
(1958), Danard (1964; 1966), Dutton and Johnson (1967), Palmen 
and Newton (1969) , Johnson (1970) , and Bullock and Johnson (1971) 
has shown that diabatic effects, especially release of latent 
heat, are indeed important and can be sufficient to offset a 
major portion of the frictional dissipation of the storm. In 
fact, it appears that during certain stages of - its life cycle, the 
storm is self-sustaining because of the energy generated by 
diabatic heating. 

From the preceeding discussion, it is evident that much 
remains to be learned about the energy-cyclogenesis process. 

While many investigators have attacked various aspects of the 
problem, a complete description of the energy processes that are 
active during the various stages of cyclone development is lacking 
in the meteorological literature. The present research was under- 
taken in an attempt to fill this void. 

c. Objectives 

The objective of the present research was to examine the 
energy budget related to cyclogenesis in the eastern United States, 
and to establish the relative importance of the various energy 
processes. In order to accomplish this objective it was necessary 
to: 

1) Develop energy budget equations suitable for a study of 
this nature. 
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2) Develop techniques for evaluation of the various terms 
in the energy budget equations using actual data. 

3) Develop a means to display the results of the numerical 
computations so that energy processes could be related to cyclo- 
genesis. 

The details of how these steps were carried out will be 
discussed fully in the remainder of this paper. At this point it 
suffices to say that primary consideration was given to developing 
equations that employ a minimum number of assumptions and that 
simulate as realistically as possible conditions that occur in the 
atmosphere. The energy budget equations that were used in the 
present research were patterned after the ones suggested by Smith 
(1970) . Also, the computational scheme used to evaluate the terms 
in energy budget equations was designed so that it would incor- 
porate all of the resolution contained in the rawinsonde 
observations over the eastern United States. 
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2. THEORETICAL DEVELOPMENT 


When the energy equations are applied to global atmospheric 
processes, the resulting expressions are somewhat simplified 
because there are no boundary integrals. This is because an 
integral of the form X i ' j " vanishes when the volume 
over which the integration takes place encloses the entire 
atmosphere. 

As shown by Wiin-Nielsen (1968), the following relation holds 
in this situation: 

f * f sL. ~ f ^IrodS/t 

dfc J v * dt ^ ^ -V»- Jit ' 

Here b is any arbitrary scalar defined in the atmosphere, £ is 
the density, and V* is the volume that encloses the entire atmo- 
sphere. If, however, one considers a volume (V) extending over 
a limited region of the atmosphere, (1) no longer holds and the 
following expression must be used: 

=Jp||dV 4 ’ (2) 

In the above expression V is the three-dimensional wind vector, 

V is the three-dimensional del operator, and V is assumed to be 
independent of time. 
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The internal, gravitational potential, and kinetic energies 
for an atmospheric volume, V, are given by 

X = C C^TdV", (3) 

Jy 


E d\J * ? (4) 

and 

K <5> 

where is the density, C v is specific heat at constant volume, 
_2 

g = 980 cm sec , and T is the virtual temperature. 

When (2) is applied to (3), the result for the time rate-of- 
change of internal energy is 

= fc.u.pdTolV" - fv-(fC lr Tv)av" («> 

An expression for C v dT/dt can be obtained by combining the 
thermodynamic equation with the continuity equation to produce 



Cu- dT - H -JLV'V, 
dt it ? 


(7) 
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where dH/dt is the diabatic heating per unit mass, and p is pres- 
sure. Upon substitution of (7) into (6), we have 



In a similar manner when (2) is applied to (4), the following 
equation results. 



where w = dz/dt is the vertical component of V. 
Finally, when (2) is applied to (5), we have 


|^ m i?tCV) av ’-J7'( f, S yv ) av » (10) 

which may be rewritten by substitution of the vector equation of 
motion to produce, 

_J vvM v '-X ? t or<iv '' 

+ J V-TF<*V' - j <U) 

where IF is the three-dimensional frictional force. 

As a matter of interest, the corresponding expressions for the 
entire atmosphere, as developed by Wiin-Nielsen (1968) are. 
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H 


(12) 


<ai = r p<^M a \ j " — r ^V'Vdiv 

*“'y * v ^ 




dP = f ^ur4V" 5 < 13 > 

Jy« 

and 

cAK a — f V # ^ d V — f +■ Tv • IF <AV • 

ot£ ^v* A/* *v* 

These expressions are obviously simpler than (8), (9), and (11) 
because of the absence of boundary integrals. 

Equations (8), (9), and (11) are the ones suggested by Smith 
(1970) and were the basis for the present research. The advantage 
of these equations is that they are derived from the thermodynamic 
equation and the equations of motion, and the individual terms 
should represent the basic physical processes that cause changes 
in I, P, and K within a specified closed volume. Before they are 
used in practice, however, it is necessary to make some modifi- 
cations. 

The first modification that must be made is to apply the 
hydrostatic approximation. This is necessary simply because the 
hydrostatic approximation is already incorporated in the rawinsonde 
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observations that are used to evaluate the various terms in (8), 
, (9), and (11). 

When the hydrostatic approximation is invoked, the vertical 
component of K, J PSjf'cW " ^ becomes invariant, and the 
expressions corresponding to (8), (9), and (11) become 


_ jj^ur R dV - - jv-^Cu-TV^V* 


(15) 




(16) 


and 


3Kh 


= - jv-(fVK V )<w’ <l7) 


where 




(18) 


is the horizontal component of K. 

Suitable justification for neglecting the vertical portion 
of K has been presented by Wiin-Nielsen (1968) and Smith (1970). 
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As pointed out by Richardson (1922), the vertical velocity is no 
longer a free variable (i.e., one that cannot be derived from 
basic parameters) in a hydrostatic atmosphere. The symbol w^ is 
used instead of w for vertical velocity in (15) and (16) to keep 
this distinction clear. 

The differences between a hydrostatic atmosphere and a non- 
hydrostatic atmosphere become clear when equations (12), (13), 
and (14) are written in symbolic form following the notation of 
Wiin-Nielsen (1968). First, for the non- hydrostatic case. 


dr > g-ci} - C(i,kV 
Jp ’ 

dP = C(K.P), 

dp 

and 

4K = C(x,k') -C(K,P) -D(k). 

dp ’ 

and, for the hydrostatic case, 

dr = g-Ci) - C(i,K h ) - C(r.p), 
dP 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 
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( 23 ) 


dip = C(x,P) 

at 

and 

i^Kk = Cfr.Kh 1 ) - d(Kk) . w 

At 

In (19) to (24), G(I) refers to generation of internal energy by 
diabatic processes, C(A,B) refers to a conversion of energy from 
form A to form B, and D is dissipation. In (20) and (21), the 
term C(K, P) appears with opposite signs, which means that in a 
non-hydrostatic atmosphere the kinetic energy may be influenced 
directly by changes in potential energy. On the other hand, this 
direct relationship between kinetic and potential energy does not 
exist in a hydrostatic atmosphere, as can be seen by an examination 
of (23) and (24). This is because in the hydrostatic atmosphere 
the potential and internal energies must remain in a constant ratio, 
as shown by Petterssen (1956); i.e., 

P * £L X . < 25 > 

Cu* 

Although the atmosphere is never in exact hydrostatic balance, 
the energetics of the real atmosphere are effectively equivalent to 
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the energetics of a hydrostatic atmosphere. This is true because 
the time required for hydrostatic adjustment is small compared to 
the time scales of the large organized motions (Johnson, 1970). 

The interpretation of the terms in (19) to (24) has been 
given by Wiin-Nielsen (1968). However, Lettau (1954) shows that 
simple mathematical manipulation can produce several different 
conversion terms in an open system. Note that the appearance of 
the terms C(K, P) in (20) and (21) with opposite signs, a necessary 
condition for an energy conversion process, is possible only 
because, for a volume enclosing the entire atmosphere, 

C t>v-volv" = -fv-vfidlv? <26) 

V -V 

However, for a limited volume, the following equation holds, 


f +• jV-V^dv" * (V- ^VdV", 

s “'v 'S\J jy 


(27) 


Smith, citing work by Pfeffer (1957) states that unless the 
boundary transport term, J S 7, | =>VdV" , is identically 
zero, neither of the two terms on the left-hand side of (27) can 
be regarded as specifying the rate of conversion between I and K. 
Smith overcomes this difficulty by writing 


fv^VdV" = 0 = (V-|3Vdv"+ fv-^VdV" 
X,<* Jv Jw' 


(28) 
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where V* and V are as defined previously, and V' is the atmospheric 
volume surrounding V, i.e.. 


V 


/ 



(29) 


Thus, the boundary work within V is balanced outside of V by that 
done within V*. This permits the interpretation of the term, 

, in (27) as a transformation of internal energy 

Jv 

within V to internal energy in the surrounding volume, V'. The 
obvious advantage of this manipulation is that it permits the 
study of the interaction of V and the surrounding atmosphere. 

When (27) is substituted into (8) along with the hydrostatic 
approximation, the following expression results: 



dlV * — (V • dV * — (V* (oCvTV)<lV* < 3 °) 

J v -V 

Finally, by defining TT = P + I, we can combine (9) and (30) to 
form 



vf>vdv" - Jv- ( P c t> -rv)dv'' (31) 

where Cp is the specific heat at constant pressure. 
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For computational purposes it is convenient to express (3), 
(4), and (5) in an x,y, p, t-system, since pressure usually is used 
as the vertical coordinate rather than height. 

We can write the total potential energy (TT*), by assuming 
hydrostatic equilibrium and using (25) as 

b « 

it * £l r r * • < 32 > 

t J A K 

In a similar manner, the horizontal component of K becomes 


K h * f t> ' V 2 -V t 
7 2 

where A is the horizontal cross-sectional area of V, A * is a dummy 
variable of integration corresponding to A, and p-^ and p2 are the 
pressures at the bottom and top of an unspecified layer in the 
atmosphere . 

Equations (17) and (31) can be transformed to an x,y,p,t 
coordinate system through the use of the following relationships 
(Thompson, 1961) , 


(HX ■ * If (1*1 ; > 

(§§>. • (!$)► - » 
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(34c) 


a<j> = d<}> ab 

B2 at> 91 

and 


=fV ^ 3 


(34d) 


where 4> is any scalar quantity. Upon substitution of (34) into 
(17) and (31) , the final energy budget equations become, 

+|j > Cr-w)]dH A * 

- f f ‘[ V l- • T ' V * + S ( T£ °)] d H«' } (35) 

l J A J p v r 
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and 




t* 


+ tU'> v *- r ^ dA ” 

< 36> 

Pj *• 

where (O = dp/dt. 

The interpretation of the terms in (35) and (36) is as 
follows. The first term on the right-hand side of (35) represents 
changes in total potential energy caused by diabatic processes. 

The second term on the right-hand side of (35) represents the rate 
at which total potential energy is changed due to cross-contour or 
non-geostrophic flow. This term also appears on the right-hand 
side of (36), but with an opposite sign, and thus may be regarded 
as the term that represents conversion between total potential and 
kinetic energy. The third term on the right-hand side of (35) cor- 
responds to the first term in (29), and as shown previously this 
term represents the conversion of total potential energy within V 
to total potential energy in the surrounding atmosphere. The last 
terms in (35) and (36) represent the flux of total potential and 
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kinetic energy through the boundaries of V. Finally, the second 
-term on the right-hand side of (36) represents changes in kinetic 
energy due to frictional forces. 


In the remainder of the text, the following symbolic notation 
will be used when referring to the terms described above: 


G(*rr) 
c.(«r .y 
c«r,Tt') 

H(ir) 

VF(TT) 

Hey 

vF(y 

D(K h ) 


Changes in total potential energy due to 
diabatic processes. 

Conversion of total potential energy into 
kinetic energy. 

Conversion of total potential energy within 
a volume to potential energy in the surrounding 
volume. 

Horizontal component of the boundary flux of 
total potential energy. 

Vertical component of the boundary flux of 
total potential energy. 

Horizontal component of the boundary flux of 
kinetic energy. 

Vertical component of the boundary flux of 
kinetic energy. 

Dissipation of kinetic energy. 


Equations (35) and (36) are the energy budget equations that 
were used to investigate the relationships between energy processes 
and cyclogenesis. The various terms in (35) and (36) were evaluated 
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with actual rawinsonde observations in a manner that will be 
explained in the next section. 
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3 . ANALYTICAL PROCEDURES 
a. The grid system 

Computation of the various terms in the energy budget 
equations is facilitated if the data points are uniformly distri- 
buted in space. Since the rawinsonde observations used in this 
study are arranged in a random manner, it was desirable to inter- 
polate the observations to equally spaced grid points. The grid 
system used in this study is shown in Fig. 1. 



Fig. 1. Grid used for numerical computations. 


25 



The grid spacing is \ in. on a 1:12,500,000 conformal conic pro- 
jection with true latitudes at 30 and 60 deg N. The conformal 
conic projection was selected because the map scale factor varies 
only slightly (0.98 to 1,02) over the area of interest. This 
slight variation was neglected, and a value of 1.0 was used in this 
study. The %-in. grid spacing is convenient also because it 
enables one to produce printed output from a computer that is the 
same scale as the 1:12,500,000 conformal conic projection. 

The %-in. spacing of grid points corresponds to a distance of 
approximately 158 km between grid points. This distance is more 
than adequate to resolve any of the information contained in the 
rawinsonde observations. In fact, in a very recent study by Barr 
et al. (1971), it is shown theoretically that a 169-km grid inter- 
val incorporates as much detail as can be justified from the 
existing rawinsonde network over the United States. 

b. Interpolation procedures 

The rawinsonde and surface observations were interpolated to 
the grid using an objective analysis scheme developed by Barnes 
(1964). This procedure is commonly referred to as successive cor- 
rections to a first-guess field. First-guess fields usually are 
derived from forecast fields, persistence, climatology, or various 
combinations of the three. Since none of these was available, zero 
was used as the first guess in this study. This has a certain 
disadvantage because forecast fields can be used to detect obvious 
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errors in the data. Since forecast fields were not available, ex- 
tensive hand checking of the data was required to eliminate errors. 

The steps in the analysis procedure are as follows. The 

th 

observations are scanned in order, and for the k observation, 

the error (eS between the observation and the first-guess analysis 

at the observation location is computed. The value of the analysis 

at the observation location is found by bi-linear interpolation. 

This error then influences a number of grid points in the vicinity 

of the observation location with the amount of influence or weight 

varying with the distance between the grid point and the observa- 

k 

tion location. The weight (w„) for grid point (i,j) is given by 


-4-r 2 / 


us 




= e 


where r is the distance from the k^ observation to the grid point, 
and SR is the scan radius. A scan radius of 4 grid units was used 
in this study. Barnes has shown this value to be optimum for the 
rawinsonde density of the eastern United States. Only the grid 
points within a distance SR from the observation point are influ- 
emced. w„ has a value of unity when the grid and observation 
points coincide, and approaches zero as the distance from the 
observation point increases. 

Any particular grid point may be influenced by several obser- 
vations, and after all of the observations have been considered, 
the total correction (c„) for each grid point is determined by 
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where N is the number of observations influencing grid point (i,j). 
The total, corrections are added to the first-guess field to provide 
a new first-guess field. 

The above procedure is repeated until there is close agreement 
between the observations and the analysis. Barnes has demonstrated 
mathematically the convergence of this procedure, and shows that 
the number of iterations for convergence depends on the amount of 
information contained in the observations. Four iterations were 
found to be sufficient for the purposes of this study. 

The choice of this technique over the more commonly used 
interpolation procedure of Cressman (1959) was motivated by two 
factors. The Cressman technique makes use of a decreasing scan 
radius on each successive iteration, whereas the scan radius is 
held fixed in the Barnes technique. It can be shown that de- 
creasing the scan radius will result in an amplification of the 
errors in the observations with each successive iteration. The 
fixed scan radius eliminates this problem and the errors in analy- 
sis due to errors in observations are constant from one iteration 
to the next. If the Cressman technique is used, extensive 
smoothing must be employed to eliminate noise introduced by the 
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errors in the observations. The smoothing in turn destroys some 
of the circulation centers. 

The Barnes technique also provides a superior representation 
of gradients. The main advantage of the Cressman technique is that 
it requires fewer computations because, after the first iteration, 
fewer grid points are influenced by each observation than with the 
Barnes technique. 

c. Description of analyzed fields 

Gridded analyses of height, temperature, and wind components 
were made for the 1000-, 850-, 700-, 500- , 400-, 300-, 200-, and 
100-mb levels. The analyses then were interpolated vertically to 
produce ten levels of information at 1000, 900, 800, 700, 600, 500, 
400, 300, 200, and 100 mb. Cubic interpolation was used to arrive 
at the temperature and wind components at the interpolated levels 
of 900, 800, and 600 mb. The heights at these levels were computed 
using the hypsometric equation. 

The data used to produce the analyses consisted of 0000 and 
1200 GMT rawinsonde and surface observations. Whenever possible, 
observations outside of the grid were used so that the boundary grid 
points would be just as reliable as the interior points. 

The analyzed fields were smoothed using a nine-point filter 
with a smoothing weight of 0.5 as described by Shuman (1957). This 
smoothing filter is defined by 


29 


A. * Ac + i V(|-V)(A, •‘A j +A, + A+-4A 0 ) 

"*■ ( As + A t + A*j + ^ A„) 

where A represents the field value at the various points indicated 
by the subscripts shown in Fig. 2, and V is the smoothing weight. 

Barr et al. (1971) have shown 
that wave lengths of less than ap- 
proximately 1700 km cannot be 
represented by the data from the 
existing rawinsonde network over 
the United States. The filter 
described above retains virtually 
1007o of the amplitude of wave lengths greater than 1500 km while 
eliminating the smaller wave lengths that are due to spurious small- 
scale fluctuations in the observations. 

d. Computation of vertical velocities 

The choice of a method for the computation of vertical velo- 
cities was not easy because none of the known methods yields 
entirely satisfactory results. The kinematic method was used in 
this study with the assumption that 60 = 0 at p = 1000 mb. This 
method was chosen because it uses actual winds and involves the 
least stringent assumptions. 


6 2 5 

3 0 1 

7 4 8 

Fig. 2. Nine-point 
smoothing mesh. 
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Two other commonly-used methods that were considered are the 
'adiabatic method and the G* -equation method. The adiabatic method 
was discarded because it produces erroneous results when applied 
to energy conversions due to the omission of diabatic processes 
(Wiin- Nielsen, 1964). Solution of the M-equation requires that CO 
be specified on the boundaries of the grid, and these boundary 
values can affect significantly the interior values because the 
grid used in this study covers a relatively small area. Other 
methods that use the quasi-geostrophic assumption were not consid- 
ered because there is a significant difference between energy 
transformations in an atmosphere that is in geostrophic balance 
versus one in ageostrophic balance. 

e. Determination of the dissipation of kinetic energy 

The frictional dissipation in the boundary layer was computed 
through the use a relationship given by Lettau (1959): 

D(Kk) ~ Cost (37) 

where ^4 

C = -fe.' [jL "Uj. - I.8fc5l , 


and Ug is the geostrophic wind speed at 1000 mb, is the angu- 
lar departure from the geostrophic wind and shear stress vectors. 


31 



/ 

is the von Karman constant, f is the Coriolis parameter, and 
Zq is the roughness parameter. The values for % , k, and Zq were 
taken to be 20 deg, 0.4, and 50 cm, respectively. These values 
have been suggested as representative for the eastern United States 
by Bullock and Johnson (1971). 

The dissipation of kinetic energy in the free atmosphere has 
been computed by Holopainen (1963), Kung (1966), and others as the 
residual term in equations similar to (36). However, these compu- 
tations were averaged over several months, and no estimates of 
free atmospheric dissipation in individual storms are available. 

As the results of this study will indicate, it is doubtful whether 
the various terms in (36) can be computed with sufficient accuracy 
for a single storm to determine dissipation as the residual term. 

An approximate expression for the dissipation term in the 
free atmosphere has been presented by Palmen and Newton (1969) as 

D(K^ = { 

where /*>> is a coefficient of lateral mixing. However, is 

not known with sufficient accuracy to make the above expression 
usable. 


32 


f . Numerical evaluation of integrals 


The integrands on the right-hand side of (32), (33), (35), and 
(36) were evaluated at each of the levels between 1000 and 100 mb. 
In general, centered differences were used except for the evalua- 
tion of time derivatives where 12-hour backward differences were 
used. The use of centered differences reduced the original 18 x 18 
grid to 16 x 16. 

Vertical integrations then were carried out through the use 
of the trapezoidal rule for the layers 1000-700 mb, 700-400 mb, and 
400-100 mb. The partial integrals for the three layers were summed 
to obtain values for the entire depth 1000-100 mb. 

O 

Units of joules m were used to express kinetic and total 
potential energy. Units of watts m were used for terms involving 
conversions of energy. Numerical computations were performed on 
the IBM 360/65 computer at Texas A&M University. 


33 


4. DESCRIPTION OF SYNOPTIC CASES 


a. Selection criteria 

Examples of cyclone development over the eastern United States 
are not difficult to find since cyclogenesis is a common occurrence 
in this area, particularly during the winter months. The primary 
consideration for the selection of synoptic cases for study was 
that the various stages of cyclogenesis had to occur while the 
cyclone was within the confines of the grid shown in Fig. 1 (p. 25). 
Only cases that occurred within the last 5 years were considered 
because of availability of data at Texas A&M University. 

The two synoptic cases that were chosen now will be discussed 
briefly. These cases have many similarities, but as will be 
shown later, the energy processes and the mechanisms responsible 
for cyclogenesis in the two cases were quite different. 

b. Case I: 1-4 November 1966 

The first indication of cyclogenesis at the surface occurred 
at 1200 GMT on 1 November 1966 as a weak wave developed on the 
frontal system through the Gulf of Mexico (Fig. 3a). The low- 
pressure center associated with this wave moved northeastward into 
Florida, and by 0000 GMT on 2 November 1966, a closed circulation 
was evident at the surface (Fig. 4a). The low center then moved 
in a north- to-northwestward direction toward the Great Lakes, and 
rapid deepening occurred until near 1200 GMT on 3 November 1966 
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(Fig. 6a) at which time the low reached its greatest intensity 
(see Figs. 5a, 6a, and 7a). By 0000 GMT on 4 November 1966, the 
low center had moved north of the Great Lakes, and was dissipating 
slowly (Fig. 8a). 

Upper- air charts for selected levels also are shown in Figs. 

3 to 8 to present a three-dimensional picture of the atmosphere 
during the period of cyclogenesis. An intense trough can be noted 
to the west of the surface low at all levels above the surface 
during the period when cyclogenesis was occurring. The southern 
portion of the trough moved eastward and sharpened while the 
northern portion remained relatively stationary. 

The approximate location of the jet stream (the position of 
the 200-mb maximum winds) is shown on the surface charts in Figs. 

3 to 8. Wind speeds in excess of 70 m sec ^ can be noted along 

the east coast of the United States in Fig. 6a. The implica- 
tions of this jet stream will be discussed later, and evidence will 
be presented that indicates that this jet stream provided the 
initial impulse of kinetic energy to the cyclone. 

This synoptic case is somewhat unusual because the track of 
the cyclone is farther west than usual. Normally, in a synoptic 
situation of this type, intense cyclogenesis will occur off the 
coast of the Carolinas, and the low center will move northward, 
remaining just off the east coast of the United States. 
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c. 500-nib chart. 


d. 200 -mb chart. 


Fig. 6. Synoptic maps for 0000 GMT, 3 November 1966. (Dashed 
lines on surface map represent 500-mb vorticity (10"^ sec"^-) . 

Dashed lines on upper- level maps are values of CD (10 - ^ mb see - •*•).) 
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a. Surface chart. 


b. 800-mb chart. 
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c. 500-mb chart. 


d. 200-mb chart. 


Fig. 8. Synoptic maps for 0000 GMT, 4 November 1966. (Dashed 
lines on surface map represent 500-mb vorticity (10"^ sec - ^). 

Dashed lines on upper- level maps are values of 03 (10~^ mb sec"l).) 





c. Case II: 5-8 December 1969 

At the beginning of this period, a weak low-pressure center 
was evident at the surface in eastern New Mexico (Fig. 9a). The 
low then moved southeastward into Texas (Figs. 10a and 11a) with 
little change in intensity. The low center then turned and moved 
northward (Figs. 12a to 14a) and the circulation around the center 
intensified considerably. By 1200 GMT, 8 December 1969, the low 
center had reached Lake Superior, and had entered the dissipation 
stage (Fig. 15a). The remarkable feature of this low is that the 
central pressure remained steady at approximately 998 mb throughout 
the period. 

The associated upper-air charts are shown in Figs. 9 to 15. 

As the low intensified at the surface, the upper-level trough also 
intensified, and by the end of the period a closed circulation can 
be noted at all levels (Figs. 14 and 15). This can be contrasted 
with the first case where this closed circulation existed only in 
the lower troposphere even when the low had reached its maximum 
development. 

Another interesting feature of this case is the organized 
patterns of vertical motion that can be seen in Figs. 11 to 13. 
These patterns were not nearly as evident in the first case. 

The role of these organized patterns of vertical motion in pro- 
ducing energy transformations and contributing to cyclogenesis 
will be discussed in detail in the next section. 
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a. Surface chart. 


b. 800 -mb chart. 
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d. 200-mb chart. 


Fig. 9* Synoptic maps for 1200 GMT, 5 December 1969. (Dashed 
lines on surface map represent 500-mb vorticity (10“5 sec~l). 

Dashed lines on upper- level maps are values of CO (10~^ mb sec"-^),) 
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5. DISCUSSION OF RESULTS 


a. Kinetic energy budget excluding dissipation 

The terms of the kinetic energy budget equation (Eq. 36), 
except for dissipation, are shown in Figs. 16 to 20 for Case I 
and in Figs. 21 to 27 for Case II. The solid lines represent 
the conversion of potential energy into kinetic energy, and 
the dashed lines represent the horizontal component of the 
boundary flux term, H(K^). The stippled regions denote areas 

where the vertical component of the boundary flux term is greater 

-2 

than 20 watts m . The position of the low center and its asso- 
ciated frontal system is shown in the upper left frame in each of 
these figures. In the discussion that follows, layer 1 will be 
used to refer to the 1000- to 700-mb layer, layer 2 to the 700- 
to 400-mb layer, and layer 3 to the 400- to 100-mb layer. 

1) Conversion terms - Case I. - At 1200 GMT on 1 November 
1966, the significant item to note is the large conversion values 
in layers 2 or 3 over the southeastern United States (Fig. 16b 
and 16c). The maximum conversion rates in layer 2 appear just 
to the west of those in layer 3, atypical of the vertical 
tilt characteristic of extratropical cyclones in the eastern 
United States. The conversion term in layer 1 appears to be 
insignificant at this time. 

In Fig. 17, which is for 0000 GMT on 2 November 1966, the 
large conversion rates still are apparent, but the patterns have 
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c. 400-100 mb. 


Fig. 17. Synoptic maps of kinetic energy budget terms for 0000 
GMT, 2 November 1966. (Solid lines represent 0(^,1^); dashed lines 
represent H(Ku); and stippled areas represent regions where VF(Kft) 
is greater than 20. All values represent integrated values through 
the given layer expressed in watts m"^.) 
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a. 1000-700 mb. 


b. 700-400 mb. 



c. 400-100 mb. 


Fig. 22. Synoptic maps of kinetic energy budget terms for 0000 
GMT, 6 December 1969. (Solid lines represent CCfTSK^); dashed lines 
represent H(%) ; and stippled areas represent regions where VF(Kh) 
is greater than 20. All values represent integrated values through 
the given layer expressed in watts nf^.) 
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c. 400-100 mb. 


Fig. 23. Synoptic maps of kinetic energy budget terms for 1200 
GMT, 6 December 1969. (Solid lines represent C(Tt,K^); dashed lines 
represent H(I^) ; and stippled areas represent regions where VF(Kh) 
is greater than 20. All values represent integrated values through 
the given layer expressed in watts m -2 .) 


a. 1000-700 mb. 
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c. 400-100 mb. 


Fig. 24. Synoptic maps of kinetic energy budget terms for 0000 
GMT, 7 December 1969. (Solid lines represent CC^pjK^); dashed lines 
represent H(Kh); and stippled areas represent regions where VF(Kh) 
is greater than 20. All values represent integrated values through 
the given layer expressed in watts m“ .) 



c. 400-100 mb. 


Fig. 25. Synoptic maps of kinetic energy budget terms for 1200 
GMT, 7 December 1969. (Solid lines represent C( / JT',Kj 1 ); dashed lines 
represent H(Kh) ; and stippled areas represent regions where VF (Kh) 
is greater than 20. All values represent integrated values through 
the given layer expressed in watts m .) 
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a. 1000-700 mb. b. 700-400 mb. 



c. 400-100 mb. 


Fig. 26. Synoptic maps of kinetic energy budget terms for 0000 


GMT, 8 December 1969. (Solid lines represent CCftjK^); dashed lines 
represent H(Kh); and stippled areas represent regions where VF(Kt,) 


is greater than 20. All values represent integrated values through 
the given layer expressed in watts m" .) 





a. 1000-700 mb. 


b. 700-400 mb. 



c. 400-100 mb. 


Fig. 27. Synoptic maps of kinetic energy budget terms for 1200 
GMT, 8 December 1969. (Solid lines represent CCTTjK^); dashed lines 
represent H(Kh) ; and stippled areas represent regions where VF (Kh) 
is greater than 20. All values represent integrated values through 
the given layer expressed in watts m - ^.) 



changed somewhat. The center of maximum conversion in layer 2 
has moved northward into Indiana, and the center in layer 3 has 
split into two parts, one off the east coast of Florida, and one 
in Indiana. Conversion rates in layer 1 still are small. 

By 1200 GMT on 2 November (Fig. 18, p. 53), the conversion 

rates have increased in layer 1 to the south of the surface low 

center. The centers of positive conversion in layers 2 and 3 

have moved slowly northward and the maximum values have increased 

-2 

by about 20 watts m . The two centers in layer 3 are becoming 
separated by large negative values. 

In Fig. 19, (p. 54), which is for 0000 GMT on 3 November, 
dramatic changes have occurred in all three layers. Significant 
conversion of potential energy into kinetic energy is very promi- 
nent in layer 1 in the vicinity of the low center. In layers 2 
and 3 the large conversion rates no longer appear above the surface 
low. The only large conversion values now appear over Florida in 
association with the jet stream. 

In the final figure in this series. Fig. 20 (p. 55), the 
low has reached its maximum intensity and has begun to dissipate 
slowly. The conversion term has become negative in the vicinity 
of the low center thus indicating a conversion of kinetic energy 
into potential energy in layer 1 of the cyclone. In layers 2 and 
3 the conversion term is small above the low center. 
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2) Boundary flux terms - Case I. - At 1200 GMT on 1 
November, the horizontal component of the boundary flux term is 
large and negative in layer 3 over the southeastern United States 
thereby indicating that most of the kinetic energy that is being 
generated in this region is being transported out of the region. 
However, in the extreme southeastern United States, the vertical 
component of the boundary flux is large and positive which implies 
that kinetic energy is being concentrated in layer 3. Since the 
vertical motion is downward in this area (see Fig. 3d, p. 35 ), 
there is also a considerable downward flux of kinetic energy in 
northern Florida and southern Alabama during this time. It is 
significant to note that this is the area in which cyclogenesis 
will begin to occur in about 12 hours. The horizontal component 
of the boundary flux term is large and negative in layer 2 in the 
east central United States, matching exactly the positive conver- 
sion rates present in that area (Fig. 16b, p. 51). The vertical 
component of the boundary flux term is positive in layer 2 

throughout the southeastern United States, but since the values 

-2 

are all less than 20 watts m , this area does not appear in 
Fig. 16b (p. 51). Finally, both components of the boundary flux 
term appear to be insignificant in layer 1 . 

By 1200 GMT on 2 November (Fig. 18, p. 53), the horizontal 
component of the boundary flux term in layer 3 has increased over 
the southeastern United States to the point where more than 
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five times as much kinetic energy is being removed horizontally 

as is being generated. The vertical component of the boundary 

-2 

flux term in layer 2 has increased to more than 60 watts m in 
the southeastern United States. Coupled with the downward verti- 
cal motion that is evident in Fig. 5c (p. 38), it is apparent that 
there is a significant downward flux of kinetic energy through 
layer 2 in this area. This represents a lag of' about 24 hours 
from the time that the large downward flux was evident in layer 3. 

In the final two figures in this series. Figs. 19 and 20 
(pp. 54 and 55), a general decrease in the values of the boundary 
flux term can be noted in layers 2 and 3, except in the vicinity 
of the jet stream where both the horizontal and vertical compo- 
nents of the boundary flux term remain large. In Fig. 19 (p. 54), 
the vertical component of the boundary flux term is large in 
layers 2 and 3 over Ohio. From Fig. 6 (p. 39) it can be seen that 
the vertical motion is upward in this area, thereby implying that 
kinetic energy in this region is being transported upward now 
rather than downward as on 1 and 2 November. These areas still 
are apparent in Fig. 20 (p. 55), but have shifted to the north 
and east and have become smaller in size. The horizontal component 
of the boundary flux term in layer 1 remained relatively insig- 
nificant until 1200 GMT on 3 November. At this time it can be 
seen in Fig. 20 (p. 55) that kinetic energy is being concentrated 
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by this term in layer 1 on the western side of the low, while 

the reverse is true on the eastern side of the low. 

3) Summary of kinetic energy budget for Case I . - At this 
point it would be appropriate to reconstruct the steps that occur- 
red during the development of this cyclone. Prior to the initia- 
tion of cyclogenesis at the surface, kinetic energy was being 
generated in the upper levels of the atmosphere over northern 
Florida. A significant portion of this energy was transported down- 
ward through the atmosphere where it supplied kinetic energy to the 
lower levels. Note the rapid increase in the 500-mb vorticity values 
in the Gulf of Mexico from 1200 GMT on 1 November 1966 (Fig. 3a, p.35) 
to 0000 GMT on 2 November 1966 (Fig. 4a, p. 36) which indicates 

the strengthening of the 500-mb circulation just prior to cyclo- 
genesis at the surface. As the cyclone began to develop, vertical 
velocities were generated that provided the direct circulation neces- 
sary for the conversion of potential energy into kinetic energy in 
the low levels of the cyclone. As the cyclone reached the mature 
stage, the upward vertical motions extended into the cold air and 
resulted in a conversion of kinetic energy into potential energy. 

The modification of the cyclone by frictional dissipation and 
diabatic processes will be discussed later. 

4) Conversion term - Case II . - The kinetic energy budget 

for the second synoptic case is shown starting with Fig. 21 (p. 56) 
which presents conditions for 1200 GMT on 5 December 1969. The first 
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item to note in Fig. 21 (p. 56) is the presence of areas of sig- 
nificant conversion of potential energy to kinetic energy in all 
three layers to the north of the warm front through southwestern 
Texas. These positive values of conversion are nearly balanced in 
layers 2 and 3 by the negative values which appear to the south of 
the front. 

By 0000 GMT on 6 December (Fig. 22, p. 57), the conversion 
pattern in layer 1 has changed so that there is now positive con- 
version on both sides of the warm front and to the north of the low 
center. The conversion values in layer 3 have become much smaller 
over Texas, while those in layer 2 have remained unchanged. 

During the next 24 hours a significant change occurred in the 
conversion term. This change can be seen in Figs. 23 and 24 (pp, 58 
and 59) which depict conditions for 1200 GMT on 6 December and 0000 
GMT on 7 December. The important change that occurred was the 
appearance of negative conversion values in layers 1 and 2 in the 
vicinity of the cyclone center. Apparently an indirect circulation 
(cold air rising) existed in this region of the cyclone that 
resulted in a conversion of kinetic energy into potential energy. 

The area of upward vertical motion to the northeast of the cyclone 
center (Fig. 11, p. 45) seems to substantiate this idea. 

By 1200 GMT on 7 December (Fig. 25, p. 60), the low center has 
moved into southern Illinois and areas of positive conversion 
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appear on all sides of the cyclone in layer 1. In layers 2 and 
3, however, small negative conversion rates are present above 
the cyclone center. 

5) Boundary flux term - Case II . - In Fig. 21 (p . 56), the 
horizontal component of the boundary flux term is large and 
positive in layers 1, 2, and 3 over southeastern Texas and the 
Gulf of Mexico, thereby indicating an increase in the concentration 
of kinetic energy at all levels in this region. The vertical 
component of the boundary flux term also is positive over this 
region, and by referring to Fig. 9d (p . 43), one can see that the 
vertical motion is upward in this area from which it can be infer- 
red that kinetic energy is being transported upward in the high 
levels of the cyclone. 

The most noticeable change during the next 12 hours is the 
appearance of large negative values of the horizontal component 
of the boundary flux term in western Texas (Fig. 22, p. 57). The 
vertical component of the boundary flux term still is large and 
positive in layer 3 over Texas, but the center has shifted north- 
ward. In layer 2, two areas can be seen where the vertical 
component of the boundary flux term is large and negative. Since 
the vertical motion generally is upward in this area (Fig. 10c, 
p. 44), it appears that kinetic energy now is being transported 
upward through both layers 2 and 3. 
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In Figs. 23 and 24 (pp. 58 and 59), the vertical component 
of the boundary flux term is much more apparent over the middle 
and high levels of the cyclone. An extensive area of upward 
vertical motion at 800 and 500 mb to the east of the cyclone center 
can be seen in Figs. 11 and 12 (pp. 45 and 46). Upward vertical 
motion also exists at 100 mb, but the magnitude is smaller. This 
pattern of vertical motion results in an upward flux of kinetic 
energy that produces the increase in intensity of the jet stream, 
as shown in Figs. 11a and 12a (pp. 45 and 46). In both Figs. 23 
and 24 (pp. 58 and 59), the horizontal component of the boundary 
flux term is large and positive in layer 1 to the northeast of the 
cyclone center thus indicating an increase in the concentration of 
kinetic energy in this region. 

By 1200 GMT on 7 December (Fig. 25, p. 60), the boundary flux 
terms indicate horizontal convergence of kinetic energy in layers 
1 and 2 on the eastern side of the low center. There still is a 
significant vertical transport of kinetic energy from low to high 
levels of the cyclone, but the area of vertical transport had 
decreased. 

In Figs. 26 and 27 (pp. 61 and 62), a general decrease in 
both components of the boundary flux terms can be noted in all 
layers of the cyclone. The only significant values appear to the 
south of the cyclone center in association with the jet stream. 
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6) Summary of kinetic energy budget for Case II . - Signifi- 
cant low- and middle- level conversion rates were apparent throughout 
the development of the cyclone in Case II. Boundary fluxes then 
were instrumental in producing an increasing concentration of 
kinetic energy in the eastern positions of the cyclone, partic- 
ularly in the low and middle levels. Part of this kinetic energy 
was transported upward into the jet stream level a uu then trans- 
ported horizontally out of the region. As will be seen later, 

part of this kinetic energy also was dissipated by frictional 
forces. As the conversion rates decreased late in the life of the 
cyclone, the gradients of kinetic energy decreased throughout the 
cyclone and resulted in a decrease in both components of the 
boundary flux terms. 

7) Comparison of kinetic e nergy budget in Cases _! and II . - 
Case I was characterized by large positive conversion rates in the 
upper levels of the cyclone during the initial development stage. 

A significant portion of the kinetic energy generated at the high 
levels was transported downward through the atmosphere, thus 
supplying a source of energy for the lower levels. As the develop- 
ment of the cyclone continued, the upper- level conversion rates 
decreased whereas positive low-level conversion rates became more 
apparent . 

The upper- level conversion rates were much smaller in the 
second synoptic case whereas positive conversion in the low levels 
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of the cyclone could be noted throughout the development stage. 

In contrast to Case I, kinetic energy in Case II was transported 
upward through the atmosphere during the early stages of develop- 
ment; this resulted in an increase in intensity of the jet stream. 

A further comparison between the two cases is presented in Fig. 28. 
The plot of kinetic energy versus time is shown for the two 
synoptic cases in this figure. The solid lines represent the 
average kinetic energy within a 500-km radius of the cyclone center, 
while the dashed lines represent average kinetic energy within a 
1000-km radius. 

There are two features of the kinetic energy traces shown in 
Fig. 28a that should be noted. First, the solid lines show 
that the maximum kinetic energy in layer 3 occurred about 12 
hours before the maximum value was reached in layer 2. Similarly, 
the maximum kinetic energy in layer 2 occurred about 24 hours 
before the maximum value was reached in layer 1. The second sig- 
nificant feature is the differences between the solid and dashed 
lines for layers 2 and 3. In both of these layers, the maximum 
value in the dashed lines occurs about 18 hours after the one in 
the solid lines. The maximum values shown by the dashed lines 
also are larger than the corresponding ones shown by the solid 
lines . 

The above observations give additional support to the inter- 
pretation of the kinetic energy conditions in Case I that was 
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given previously. The conclusion that kinetic energy was gener- 
ated first in the high levels of the atmosphere and then transported 
downward may be inferred from the graphs in Fig. 28a (p. 72). The 
transport of kinetic energy from high to low levels of the atmosphere 
would account for the observed lag in the peaks of the kinetic 
energy traces shown in Fig. 28a (p. 72). The importance of the 
boundary flux terms in layers 2 and 3 also are demonstrated by this 
figure. That kinetic energy was transported outward from the 
center of the storm is suggested by the observed difference between 
the solid and dashed lines of Fig. 28a (p. 72). 

The corresponding kinetic energy traces for the second synoptic 
case exhibit much different characteristics. First, the differences 
between the solid and dashed lines are much less noticeable. 

Secondly, the peak in all three levels occurred at about the same 
time, 0000 GMT on 7 November 1966. The fact that the peaks in the 
three layers occurred simultaneously seems to preclude vertical 
transport of kinetic energy as being a dominant factor as in Case 
I. However, as has been pointed out, there was an upward transport 
of kinetic energy early in the period that contributed to the 
intensification of the jet stream. This increase in kinetic energy 
in layer 3 can be noted in Fig. 28b (p. 72) from 0000 GMT on 6 
December to 0000 GMT on 7 December. The similarity in the solid 
and dashed lines in Fig. 28b (p. 72) indicates that the boundary 
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flux of kinetic energy was of lesser importance in Case II than 
in Case I. 

b. Modifying influences exerted by low - level frictional 
dissipation and diabatic processes 

Synoptic maps of the diabatic heating term and the boundary 
layer dissipation term for the 1000- to 700-mb layer are presented 
in Fig. 29 for Case I, and in Fig. 30 for Case II. The solid 
lines represent values of diabatic heating computed from the 
thermodynamic equation in the form, 

+ <yv,.^ T +(<*$£ -•*) <*>. os) 

The dashed lines represent values of boundary layer dissipation, 
as computed from (37) . The location of the surface low and its 
associated frontal system also is shown to illustrate the corre- 
spondence between low-level diabatic heating and cyclone development. 

1) Case I.- At 1200 GMT on 1 November an extensive area of 
positive generation of potential energy can be noted in eastern 
Texas (Fig, 29a). This area of positive generation coincides 
with the precipitation in this area that resulted from the 
warm moist air of the Gulf of Mexico rising over the frontal 
system. The fact that the precipitation in eastern Texas corre- 
sponds to the upward vertical motion shown in Figs. 3b and 3c 
(p. 35), attests to the validity of the vertical motion 
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a. 1200 GMT, 1 November 1966 




510 15*0 \ 


b. 0000 GMT, 2 November 1966 





c. 1200 GMT, 2 November 1966 

Fig. 29. Synoptic maps of low-level frictional dissipation 
and diabatic heating for Case I. (Solid lines represent diabatic 
heating in units of 10^ watts m“ 2. Dashed lines are values of 
frictional dissipation in units of watts m - ^. Stippled areas 
denote precipitation.) 








c. 1200 GMT, 6 December 1969 



d. 0000 GMT, 7 December 1969 



Fig. 30. Synoptic maps of low-level frictional dissipation and 
diabatic heating for Case II. (Solid lines represent diabatic heating 
in units of 10^ watts m“^. Dashed lines are values of frictional 
dissipation in units of watts m“^. Strippled areas denote precipita- 
un.) 





1200 GMT, 7 December 1969 


f. 0000 GMT, 8 December 1969 



Fig. 30. (continued) 




calculations in this region. The large values of boundary layer 
dissipation in eastern Texas also should be noted. These large 
values were produced by the increased pressure gradient which 
resulted in increased wind velocity as the high-pressure ridge 
moved southward into Texas. 

In Fig. 29b (p. 75), the surface low has moved into northern 
Florida and the diabatic heating term is small over the entire 
map. There is, however, considerable dissipation in the Gulf of 
Mexico associated with the strong pressure gradient behind the 
cold front. 

By 1200 GMT on 2 November, the diabatic heating term has 

-2 

become prominent again, exceeding 600 watts m in northern Florida 

(Fig. 29c, p. 75). During this period the low was undergoing rapid 

development, and the precipitation associated with the low had 

begun to increase. The boundary layer dissipation had increased 

_2 

also to more than 20 watts m in northern Mississippi and 
Alabama. 

The fact that the area of positive generation is south of 
the low, rather than coincident with the low as would be expected, 
should be noted. The shift in the center of the area of positive 
generation away from the low center can be attributed to the manner 
in which the diabatic heating term was computed. The first term 
on the right-hand side of (38) was evaluated by subtracting tem- 
peratures which were 12 hours apart, whereas the remaining terms 
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in (37) were evaluated from data at a single time period. Since 
there was rapid movement of the low center and the cold front 
during this period, it is likely that large errors in result when 
using 12-hour backward differences. The values shown in Fig. 29c 

(p. 75) probably are more representative of average conditions 
during the 12-hour period 0000 to 1200 GMT on 2 November rather 
than at 1200 GMT on 2 November. If hourly temperature data were 
available, it would seem reasonable to expect the areas of dia- 
batic heating to coincide more closely with the low center and its 
related precipitation pattern. Indeed, in Figs. 29d and 29e 
(p. 76) this correspondence does exist as the movement of the low 
center begins to slow down. 

In Figs. 29d and 29e (p. 76) the areas of positive genera- 
tion appear just to the west and north of the low center* closely 
approximating the precipitation patterns that were present at 
these times. By 0000 GMT on 3 November, the entire western half 
of the cyclone is under the influence of large frictional dissipa- 
tion (Fig. 29d,p. 76). This is approximately the time when the 
cyclone had reached its maximum intensity. During the next 12 
hours the storm decreased in intensity and, as can be noted in 
Fig. 29e (p. 76), the area of boundary layer dissipation has de- 
creased considerably. At the same time the generation of potential 
energy by diabatic heating now exceeds 800 watts m . 
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In summary, kinetic energy was being destroyed by low-level 
frictional dissipation throughout the duration of the cyclone. The 
maximum dissipation occurred toward the rear of the storm and at 
about the time when the storm had reached its greatest intensity. 

As the cyclone reached the mature stage, a significant amount of 
potential energy was being generated by diabatic processes to the 
west and north of the cyclone center. The maximum generation 
occurred to the north of the cyclone center shortly after the 
cyclone had reached its greatest intensity. The extent to which 
the potential energy that was generated by diabatic heating was 
converted into kinetic energy cannot be ascertained precisely since 
only a small fraction of the potential energy that exists in the 
atmosphere at any particular time is ever converted into kinetic 
energy. However, the results that have been presented for Case I 
do indicate that sufficient potential energy was generated to 
offset at least a major portion of the frictional dissipation of 
the cyclone. These results also demonstrate how diabatic processes 
can serve to prolong the lifetime of an extratropical cyclone. 

2) Case II.- In contrast with Case I, the diabatic heating 
played a significant role throughout the development of the 
cyclone in Case II. The maximum contribution from diabatic pro- 
cesses appeared at 1200 GMT on 7 December (Fig. 30e, p. 78). This 
is about 12 hours after the cyclone reached its maximum intensity 
in terms of kinetic energy (see Fig. 28b, p. 72). After 1200 GMT 
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on 7 December, the diabatic heating term decreased and became 

insignificant by the end of the period (1200 GMT on 8 December 
(Fig. 30g, p. 78)). 

The frictional dissipation term was small except for the 

period 1200 GMT, 7 December to 0000 GMT, 8 December during which 

-2 

time it exceeded' 25 watts m . The areas of frictional dissipation 
consistently appear to the east of the low center in sharp contrast 
to the first case in which the frictional dissipation was concen- 
trated on the western side of the low center. 

The differences in the two cases point out the different ways 
that diabatic processes can affect cyclogenesis. In the first case 
the diabatic heating term became significant late in the development 
of the cyclone, and reached its maximum value after the cyclone had 
become fully developed. It appears that in this case diabatic 
heating served to prolong the life of the cyclone rather than con- 
tributing significantly to the cyclogenesis process. In the 
second case, however, significant diabatic heating was evident 
early in the development stage, and apparently was a contributing 
factor to the initial growth of the cyclone as well as serving to 
extend the life of the cyclone. 

c. Diabatic heating in the middle and high levels of the cyclone 
At this point only the diabatic heating for the 1000- to 
700-mb layer of the cyclone has been discussed. In Fig. 31, 
a comparison between the low- and middle- level diabatic heating 
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b. December 1969. 


Fig. 31. Diabatic heating term. (Values were obtained by 
averaging diabatic heating term within a 1000-km radius of the 
cyclone center.) 



can be seen. The low-level values are given by the dashed lines 
and the middle-level values by the solid lines. These values 
represent averages within a 1000-km radius of the cyclone center. 
In each of the synoptic cases, the low- and middle- level diabatic 
heating terms appear to be very similar. The most noticeable 
difference is that, in both synoptic cases, the diabatic heating 
is greater in the middle levels of Lae cyclone during the initial 
development. The graphs in Fig. 31 (p. 83) again point out the 
different character of the diabatic heating in the two synoptic 
cases. 

The diabatic heating term for the 400- to 100-mb layer of the 
cyclone was very erratic in both synoptic cases. The validity of 
the computations in this layer is highly questionable because of 
observational errors in the high levels of the atmosphere, partic- 
ularly at the 100-mb level. The tendency of the kinematic method 
to overestimate vertical velocities at the 100-mb level has been 
pointed out by O'Neill (1966) and O’Brien (1970). Also, the 
temperature and wind data are much less reliable at this level 
than in the lower levels of the atmosphere. The combination of 
these factors can produce large errors in the diabatic heating 
term for the 400- to 100-mb layer of the cyclone. For this reason 
it was felt that the erratic results for this layer are not 
representative of actual atmospheric conditions, and consequently 
they are not presented here in the form that they were for the 


84 


other two layers. However, the diabatic heating values for the 

400- to 100-mb layer were averaged over the life of the cyclone 

and within a 1000-km radius of the center. This procedure yielded 

_2 

values of -78 and 190 watts m for Cases I and II, respectively. 
This averaging process removes some of the spurious fluctuations 
from the data and the average values should be more representative 
of atmospheric conditions. However, there is no apparent reason 
why a negative value was obtained in Case I and a positive value 
in Case II. This difference may be simply a reflection of the 
errors in the observations. 

As a final note, diabatic heating results from release of 
latent heat during precipitation, short and long wave radiation, 
and transfer of sensible heat as the cyclone moves over heat 
sources or sinks. No attempt has been made to separate these 
components in this study. Qualitatively at least, it could be 
seen that release of latent heat provided the major contribution 
to diabatic heating in both synoptic cases because, invariably, 
large values of diabatic heating were associated with areas of 
precipitation. Estimates by Anthes and Johnson (1968) indicate 
that in the case of hurricane Hilda (1954), 777. of the generation 
of available potential energy resulted from latent heat release, 

177. by infrared cooling, and 67. by direct solar absorption. 

Sensible heat exchange at the surface of Earth was found to be a 
negligible quantity. Certain assumptions had to be made concerning 
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the amount and distribution of solar absorption and infrared 
cooling in the hurricane. The validity’ of the results hinges on 
the correctness of these assumptions. Obviously, more work is 
needed in this area. 

d. Average kinetic energy budget 

The magnitude of the terms in the kinetic energy budget 
equation are presented in Tables 1 and 2 for synoptic Cases I and 
II, respectively. The entries in these tables represent average 
values within a 1000-km radius of the cyclone center between 1000 
and 100 mb . 

In Table 1, the items of interest to note are: 

1) The large values for C(fF, K^), H(K^), and VF (K^) during 
the first three time periods. The significance of these items has 
been discussed previously. 

2) The total conversion rate of 25.3 watts m at 1200 GMT 
on 2 November compares favorably with estimates obtained by other 
investigators through the use of the &K- technique. Some difference 
should be expected because of the different expression used for 
conversion in this study. Also, the figures in Table 1 are based 
on integration to 100 mb, whereas in previous studies the integra- 
tion was terminated at 300 mb or below. 

3) In general, the boundary flux terms are of the same order 
of magnitude as the conversion term, but slightly smaller. Excep- 
tions to this can be found in layer 3 but., as already pointed out, 
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errors are more prominent in this layer and the absolute value of 
the numbers is less reliable. 


Table 1. Kinetic energy budget for synoptic Case I. Values 
in table represent averages in units of watts m~2_ within a 1000-km 
radius of the cyclone center. An explanation of the symbols used 
in the table can be found on p. 23. 


DATE 

LAYER 

c(fr,K h ) 

H <V 

VFOy 

D(K h ) 


1 

5.7 

1.4 

-0.8 

-3.3 

1200 GMT 

2 

6.6 

3.8 

-3.3 


1 Nov* 

3 

6.4 

-19.0 

7.6 



Total 

18.7 

-13.8 

3.5 



1 

4.3 

0.6 

-0.9 

-3.7 

0000 GMT 

2 

5.5 

-1.2 

3.8 


2 Nov. 

3 

14.0 

-9.0 

7.8 



Total 

23.8 

-9.6 

10.7 



1 

7.8 

1.2 

-2.1 

-6.8 

1200 GMT 

2 

5.5 

-1.5 

2.3 


2 Nov. 

3 

12.0 

-23.0 

10.7 



Total 

25.3 

-23.3 

7.3 



1 

6.5 

1.0 

-1.2 

-8.7 

0000 GMT 

2 

2.3 

-2.5 

-0.3 


3 Nov. 

3 

9.0 

-11.0 

1.5 



Total 

17.8 

-12.5 

0.0 



1 

-1.9 

1.0 

-1.4 

-6.4 

1200 GMT 

2 

2.5 

1.7 

-12.0 


3 Nov. 

3 

-2.6 

1.5 

2.9 



Total 

-2.0 

4.2 

-11.5 



1 

-1.2 

-2.1 

-0.0 

-6.2 

0000 GMT 

2 

6.3 

-2.2 

0.3 


4 Nov. 

3 

-7.4 

1.0 

-5.1 



Total 

-2.3 

-3.3 

-4.8 


In Table 

2, the 

following items should 

be noted: 


1) The 

conversion rates in 

the upper 

levels of 

the cyclone 


are much smaller than in Case I. 
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Table 2. Kinetic energy budget for synoptic Case II. Values 
in table represent averages in units of watts in”2 within a 1000-km 
radius of the cyclone center. An explanation of the symbols used 
in the table can be found on p. 23. 


DATE 

LAYER 

C(f>,K h ) 

H(K h ) 

VF(K h ) 

D(K h ) 


1 

3.8 

0.8 

-0.4 

-2.0 

1200 GMT 

2 

-0.3 

3.7 

-3.9 


5 Dec. 

3 

-0.5 

11.0 

2.4 



Total 

3.0 

15.5 

-1.9 



1 

4.9 

2.0 

-1.6 

-3.5 

0000 GMI 

O 

L. 

-4.2 

4.2 

-0.8 


6 Dec. 

3 

4.0 

-37.0 

14.0 



Total 

4.7 

-30.8 

11.6 



1 

6.4 

3.7 

-3.2 

-6.5 

1200 GMT 

2 

-4.8 

4.8 

-12.0 


6 Dec. 

3 

-0.8 

-26.0 

13.0 



Total 

0.8 

-17.5 

-2.2 



1 

7.7 

4.1 

-0.2 

-5.9 

0000 GMT 

2 

-0.4 

6.8 

-2.8 


7 Dec. 

3 

8.4 

-20.0 

4.8 



Total 

15.7 

-9.1 

1.8 



1 

8.1 

7.5 

-1.2 

-5.7 

1200 GMT 

2 

-0.8 

1.7 

0.7 


7 Dec. 

3 

4.2 

-2.3 

-1.0 



Total 

13.1 

6.7 

-1.5 



1 

2.7 

-0.7 

0.3 

-3.9 

0000 GMT 

2 

0.4 

-0.3 

-1.1 


8 Dec. 

3 

5.8 

5.1 

-2.6 



Total 

8.9 

4.1 

-3.4 



1 

3.3 

0.2 

-0.4 

-3.8 

1200 GMT 

2 

0.5 

-0.8 

-0.1 


8 Dec. 

3 

4.5 

3.4 

-0.9 



Total 

9.3 

2.8 

-1.4 


2) 

The values 

of H(K h ) for 

layer 3 probably are in error at 

0000 and 

1200 GMT on 6 December, 

and 0000 

GMT on 7 

December. The 


kinetic energy traces shown in Fig. 28b (p. 72) for these dates do 
not indicate outward transport of kinetic energy of the magnitude 
indicated by the values in Table 2. 
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3) The total conversion rate of 15,7 watts at 0000 GMT, 

7 December 1969 again compares favorably with estimates of other 
investigators. It is, however, much smaller than that for Case I; 
this attests to the fact that the development was much more spec- 
tacular in Case I than in Case II. 

It should be pointed out that average values present a some- 
what distorted picture. For example, in Case I the low-level 
conversion rates for the first four time periods are 5.7, 4.3, 7.8, 
and 6.5 watts m “2. x n Case II the conversion rates are 3.8, 4.9, 
6.4, and 7.7 watts m~2 # From these numbers it would appear that 
the two cases were very similar. However, as was shown previously, 
the magnitude of the low-level conversion rates were larger in 
Case II than in Case I even though the average values for the two 
cases are similar. This averaging process, which has been used in 
all previous studies, masks the variability of the conversion rates 
that exists in various parts of the cyclone. 

An interesting comparison can be made between this investiga- 
tion and the one by Eddy (1965). The results for Case I support 
his findings that the maximum conversion occurs early in the 
development cycle. However, the results of Case I indicate that 
the maximum conversion rates were in the upper troposphere rather 
than in the middle troposphere as Eddy found. Figure 32 shows the 
vertical profile of conversion rates for each case at the time the 
most intense cyclogenesis was occurring. 
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a. 0000 GMT, 2 November 1966 b. 0000 GMT, 7 December 1969 

Fig. 32. Vertical profiles of conversion of potential energy 
into kinetic energy. 

The profile for 0000 GMT on 2 November 1966 is very similar to 
one found by Sechrist and Dutton (1970). On the other hand, the 
profile for 0000 GMT on 7 December 1969 is more typical of findings 
by Kung (1966), particularly with respect to the sharp minimum in the 
profile in the middle layer. From Tables 1 and 2 (pp. 87 and 88) 
it is apparent also that the vertical profiles vary considerably 
with time. These examples illustrate the danger of attempting to 
generalize results based on any single case. While extratropical 
cyclones are similar in many respects, it is apparent also that each 
one has distinctive features. 
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e. Average potential energy budget 


The terms of the potential energy budget equation are pre- 
sented in Tables 3 and 4 for synoptic Cases I and II, respectively. 
The terms representing the generation of potential energy by 
diabatic processes, and the conversion of potential energy into 
kinetic energy, already have been discussed in some detail. The 
remaining terms of the potential energy budget now will be con- 
sidered. 

1) Boundary flux of potential energy . - It is apparent 
immediately in Tables 3 and 4 that the boundary flux terms are very 
large. In general, the two components, H(1T) and VF(*7T), are nearly 
equal in magnitude, but opposite in sign. Therefore, small errors 
in either component may result in a large error in the total flux 
term which is the sum of H(*fT) and VF ('IT). The magnitude of the 
flux terms also demonstrates the futility of attempting to compute 
changes in kinetic energy from corresponding changes in potential 
energy. An error of only 10% in computing the potential energy 
flux terms could result in errors in the potential energy budget 
of 700 to 800 watts m -2 , which is one to two orders of magnitude 
larger than any of the terms in the kinetic energy budget equation. 

It should be noted that the potential energy flux terms 
generally are much smaller in layer 2 than in either of the other 
two layers. This tendency can be observed in the values in Tables 
3 and 4 for both synoptic cases. This fact seems to indicate th' 
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the conditions in layer 2 were more nearly barotropic and geo- 

strophic than in the other two layers resulting in a smaller 
advection of temperature, hence a smaller advection of potential 
energy. 


Table 3. Total potential energy budget for synoptic Case I. 
Values in table represent averages in units of watts within a 
1000-km radius of the cyclone center. An explanation of the sym- 
bols used in the table can be found on p. 23. 


DATE 

LAYER 

G(1T ) 

-c(tr,K h ) 

H(-tT) 

vF(-rr) 

c(<r, it') 


1 

-120 

-5.7 

2500 

-2600 

-29 

1200 GMT 

2 

82 

-6.6 

490 

-380 

31 

1 Nov. 

3 

-350 

-6.4 

-4100 

3900 

-57 


Total 

388 

-18.7 

-1110 

920 

55 


1 

77 

-4.3 

3300 

-2500 

230 

0000 GMT 

2 

-43 

-5.5 

-580 

760 

50 

2 Nov. 

3 

-230 

-14.0 

-2000 

1800 

-57 


Total 

-196 

-23.8 

720 

60 

223 


1 

9 

-7.8 

4900 

-3600 

370 

1200 GMT 

2 

80 

-5.5 

-1200 

1700 

140 

2 Nov. 

3 

-430 

-12.0 

-2700 

2500 

-57 


Total 

-341 

-25.3 

1000 

600 

453 


1 

60 

-6.5 

1100 

-1100 

0 

0000 GMT 

2 

80 

-2.3 

-260 

400 

39 

3 Nov. 

3 

-57 

-19.0 

310 

-230 

23 


Total 

83 

- 27.8 

1150 

-930 

62 


1 

140 

1.9 

2000 

-2300 

-86 

1200 GMT 

2 

79 

-2.5 

1000 

-1500 

-143 

3 Nov. 

3 

710 

2.6 

2100 

-2000 

29 


Total 

929 

2.0 

5100 

-5800 

-200 


1 

38 

1.9 

1300 

-900 

114 

0000 GMT 

2 

-85 

-2.5 

-2600 

2000 

-172 

4 Nov. 

3 

-120 

2.6 

1000 

-1800 

-229 


Total 

-267 

2.0 

-300 

-700 

-287 
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Table 4. Total potential energy budget for synoptic Case II. 
Values in table represent averages in units of watts m - ^ within a 
1000-km radius of the cyclone center. An explanation of the sym- 


bols used 

in the table 

can be 

found on p. 23 

• 



DATE 

LAYER 

G(rr) 

-CCiT,]^) 

HC-rr) 

VF (IT) 

c (nvir' ) 


1 

-35 

-3.8 

4000 

-3100 

258 

1200 GMT 

2 

67 

0.3 

1600 

-1600 

0 

5 Dec. 

3 

300 

0.5 

1700 

-450 

358 


Total 

332 

3.0 

7300 

-4150 

616 


1 

210 

-4.9 

8600 

-7900 

200 

0000 GMT 

2 

300 

4.2 

750 

1100 

530 

6 Dec. 

3 

460 

-4.0 

-4900 

5600 

200 


Total 

970 

- 4.7 

4450 

-1200 

930 


1 

160 

-6,4 

7600 

- 7200 

114 

1200 GMT 

2 

330 

4.8 

420 

810 

352 

6 Dec. 

3 

250 

0.8 

-3200 

4500 

372 


Total 

740 

-0.8 

4820 

-1890 

838 


1 

140 

- 7.7 

3000 

-2500 

144 

0000 GMT 

2 

180 

0.4 

-370 

690 

228 

7 Dec. 

3 

-360 

-8.4 

-2100 

2900 

86 


Total 

-40 

-15.7 

530 

1090 

458 


1 

250 

-8.1 

1400 

-1100 

86 

1200 GMT 

2 

220 

0.8 

220 

-720 

-143 

7 Dec. 

3 

-1100 

-4.2 

-930 

950 

5 


Total 

-630 

-13.1 

690 

-870 

-52 


1 

130 

-2.7 

240 

290 

152 

0000 GMT 

2 

140 

-0.4 

840 

-1300 

-132 

8 Dec. 

3 

1100 

-5.8 

2400 

-1700 

200 


Total 

1370 

-8*9 

3480 

-2710 

220 


1 

-21 

-3.3 

1800 

-1400 

104 

1200 GMT 

2 

36 

-0.5 

-230 

-47 

-79 

8 Dec. 

3 

710 

-4.5 

1700 

-1800 

-29 


Total 

725 

-9.3 

3270 

-3247 

6 
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2) Interaction of the cyclone with larger scales of motion . - 
The term CdT^TT 1 ) which appears in the potential energy budget 
equation provides a unique means to study the interaction of a 
cyclone with the larger scales of motion. The average values of 
CCIT^TT 1 ) for each synoptic case are shown in Tables 3 and 4 
(pp. 92 and 93). 

In synoptic Case I, the values of C(tT,iT f ) for the entire 
layer, 1000 to 100 mb, are large and positive during the early 
periods of cyclogenesis, thus indicating that the potential energy 
of the cyclone is being increased at the expense of potential energy 
in the atmospheric volume surrounding the cyclone. After 1200 GMT 
on 2 November 1966, the values of CCtT^TT 1 ) decrease and eventually 
become negative, which indicates a reverse flow of potential energy. 

By contrast, the values of CCfTjffO f or the total layer in 
Case II are positive throughout almost the entire life of the 

cyclone. The only exception is at 1200 GMT on 7 December 1969 at 

-2 

which time a value of -52 watts m can be noted. As with Case I, 
the values do decrease with time but never reach the large negative 
values attained in Case I. The values of are much 

larger in Case II than in Case I during the early development 
stage; this suggests that the interaction of the cyclone with the 
larger scales of motion played a more important role in cyclo- 
genesis in Case II than in Case I. 
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To the author's knowledge, this is the first time that the 
term CCf]' , •fT' ’ ) has been evaluated in an energy budget study of 
this type. Therefore, the extent to which the results just 
described apply to other cyclones cannot be determined by comparison 
with other studies. 
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6. CONCLUSIONS 


A method has been presented to investigate the atmospheric 
energy budget as related to cyclogenesis. This work represents an 
extension of previous studies in that all of the terms of the 
energy budget equation were evaluated throughout the development 
period of the cyclone. Also, objectively analyzed data were used 
in the evaluation of the energy budget terms in order to minimize 
systematic errors. An objective analysis scheme is described that 
insures that all of the resolution contained in the rawinsonde 
observations is incorporated in the analyses. 

The energy budget equations used in the present research were 
patterned after those suggested by Smith (1970), but were modified 
to include the hydrostatic approximation and the use of pressure 
as the vertical coordinate rather than height. These equations 
have been shown to be advantageous because the individual terms 
represent basic physical processes which produce changes in atmo- 
spheric energy, and the equations provide a means to study the 
interaction of the cyclone with the larger scales of motion. 

Two examples of cyclogenesis over the eastern United States 
were chosen for study. One of the cases (1-4 November 1966) 
represented an example of vigorous development, while the develop- 
ment in the other case (5-8 December 1969) was more modest. While 
computational errors may be serious due to errors in the measured 
data, particularly in the higher levels of the atmosphere, it is 
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felt that the results of this study present a representative 
picture of energy processes during cyclogenesis. The major con- 
clusions that are evident from the present research are as follows: 

1) Diabatic processes can be sufficient to offset a major 
portion of the frictional dissipation of the cyclone and in some 
cases may contribute to the initial development of the cyclone. Dia- 
batic processes were significant in both of the synoptic cases that 
were considered in this paper. 

2) The downward transport of kinetic energy from the jet 
stream level can be an important source of energy for a developing 
cyclone. This factor was significant in one of the synoptic cases, 
but not in the other, which implies that downward transport of 
kinetic energy is not a necessary condition for cyclogenesis. 

3) In general, the conversion rates determined in this study 
compared favorably with those found in previous studies. However, 
it was shown that the W* -technique of computing conversion rates 
that was used in most of the previous studies may produce erroneous 
results because of the omission of boundary terms. In fact, the 
results of the present research clearly indicate that the validity 
of any energy budget study of a limited region of the atmosphere 
hinges upon the extent to which boundary fluxes of energy can be 
determined. In general, the boundary flux terms are the same order 
of magnitude as the conversion terms and cannot be neglected. 

In addition to the above conclusions, other items of interest 
to note from the present research are as follows: 
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1) The interaction of the cyclone with the surrounding atmo- 
sphere was found to be an important factor in the potential energy 
budget for both of the synoptic cases considered. 

2) When the terms of the energy budget equations are averaged 
in time and space, some of the spurious fluctuations in the data 
are eliminated, but at the same time some of the real variability 
that exists in various parts of the cyclone is masked. 

3) It is difficult, if not impossible, to draw generalized 
conclusions from the study of one or two cyclones because of the 
marked variability than can exist between energy processes in 
different cyclones. 
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7. RECOMMENDATIONS FOR FURTHER RESEARCH 


It is likely that further studies of the type presented here 
can contribute substantially to the understanding of atmospheric 
energetics as related to cyclogenesis. The approach presented in 
this study should be extended to include partitioning of the 
diabatic heating term into individual components resulting from 
latent heat release, sensible heat transfer, radiational losses, 
etc. The results concerning the interaction of cyclones with the 
larger scales of motion must be considered tentative at this point. 
Further research could determine the extent to which other cyclones 
interact with global motions. 

The main limiting factor in energy budget studies continues to 
be the lack of adequate observations, both in number and accuracy. 
Very precise measurements are required in order to determine accu- 
rately vertical motions and boundary fluxes of energy. The Global 
Atmospheric Research Program may provide a source for observations 
that are more suitable for energy budget studies than currently 
exist. These data should become available in the mid 1970' s. 
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